Observation of excited states in a graphene double quantum dot 
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We study a graphene double quantum dot in different coupling regimes. Despite the strong 
capacitive coupling between the dots, the tunnel coupling is below the experimental resolution. We 
observe additional structures inside the finite-bias triangles, part of which can be attributed to 
electronic excited dot states, while others are probably due to modulations of the transmission of 
the tunnel barriers connecting the system to source and drain leads. 
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I. INTRODUCTION 

Double quantum dot structures are promising candi- 
dates for the implementation of solid state spin qubits^^^. 
Double dots have been realized in many different mate- 
rial systems, as for example in GaAs heterostructures^, 
semiconductor nanowiresi^ ^ and carbon nanotube^^, and 
the control of individual electrons and spins has been 
achieved^^^. Graphene has been predicted to be par- 
ticularly well-suited for spin-based quantum information 
processing, because spin-orbit interaction and hyperfine 
interaction are expected to be much weaker than in the 
material systems mentioned above, leading potentially 
to much longer spin coherence timesi^ii^. Significant 
progress has been made recently in the fabrication and 
the understanding of graphene-based nanostructures, as 
for example constrictions and quantum dots^^^^. Also 
graphene double quantum dots have been demonstrated 
recently^^^^^. For example, Liu et al.— showed the pres- 
ence of excited states in a double dot created in a top- 
gated graphene nanoribbon. In this work, we demon- 
strate the presence of excited states in a side-gated 
graphene double dot structure with a different geometry, 
formed by etching the islands out of a graphene flake. We 
study in detail the coupling between the two dots for dif- 
ferent gate voltages. We show that an in-plane magnetic 
field changes the excited states spectrum. 



II. EXPERIMENTAL DETAILS 

The sample consists of a double dot structure carved 
out of a graphene flake. Mechanical exfoliation of natu- 
ral graphite flakes, followed by deposition onto a highly 
doped silicon wafer covered by 285 nm of silicon dioxide, 
is used to produce the graphene flakes. Thin flakes are 
identified with an optical microscope, and Raman spec- 
troscopy is used to make sure that the flake consists only 
of one single graphene laye r^^i^^ . The selected flake is 
contacted using electron beam lithography (EBL) and 
metal evaporation (Cr/Au). Finally it is patterned into 
the double dot structure shown in Fig. [TJa) by a second 
EBL step and reactive ion etching based on argon and 
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FIG. 1: (a) AFM image of the sample topology. The edges 
of the graphene regions are highlighted by dashed lines. The 
two dots, labelled by L and R, have a diameter of 90 nm and 
are seperated by a 30 nm wide constriction. The dots are 
connected by 20 nm wide constrictions to source and drain 
contacts (S and D). The global back gate and five in-plane 
graphene gates (CL, GL, GC, GR, CR) allow to tune the 
sample, (b) Current through the double dot measured as 
a function of Vgr and Vgl at T — 1.4 K, Vbg = Vgc = 
V, VcR = -0.1 V, VcL = -3.25 V and Vbias = -0.5 mV. 



oxygen. 

The transport experiments were carried out in a vari- 
able temperature insert at 1.4 K, and at 120 mK base 
temperature of a standard ^He/^He dilution refrigerator. 
In total, measurements of three different cool-downs are 
presented. Even if some details changed from one cool- 
down to the other, the main features presented in this 
work were present in every cool-down. 



III. RESULTS AND DISCUSSION 

Fig. mb) shows a measurement of the charge stability 
diagram, recorded at T = 1.4 K. The hexagon pattern 
characteristic for double dots is clearly visible, and uni- 
form over many double dot charge configurations. The 
current is maximal at the triple points, where the elec- 
trochemical potentials in both dots are aligned with each 
other and with the Fermi energy in the leads. These 
triple points are connected by faint lines of much smaller 
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FIG. 2: T=1.4 K, Vbias = 0.5 mV^ Ak^T (a) Current as a 
function of Vbg, all the other gates at V. (b)-(e): Current as 
a function of the voltage applied to GR and GL for different 
negative values of back gate voltage. 



current, originating from inelastic cotunneling processes. 
Along these lines, the energy level in one dot is aligned 
with the electrochemical potential in the corresponding 
lead. Such well controlled double dot behavior, some- 
times with less symmetric barriers, could be observed in 
the whole accessible range of positive back gate voltages 
(0-30 V) for all three cool-downs, as long as the barriers 
were not too closed to allow current detection. 

For negative back gate voltages, however, the situa- 
tion can be quite different. This difference can be ob- 
served in Fig. [2fa), displaying the current through the 
double dot as a function of back gate voltage at 1.4 K. 
While the current is completely suppressed for positive 
values of Vqg-, resonances can be observed at negative 
gate voltages. Charge stability diagram measurements 
give a better understanding of this region. Fig. [2fb)-(e) 
represent a selection of such measurements for four dif- 
ferent values of Vbg between -10 V and -20 V, where the 
resonances are strongest. They correspond to completely 
different situations: at Fbg = —12 V, represented in Fig. 
[2](b), the current is high along the horizontal cotunnel- 
ing lines and almost continuous across the triple points. 
Fig. [2jc) represents the opposite situation, occurring at 
Vqq = —13 V, with high current at the triple points and 
along the vertical cotunneling lines. At Vqg = — 16 V 
[Fig. [2fd)], the regions of finite current describe diago- 
nal, wavy lines, while at Vbg = — 19 V [Fig. Efe)] a clean 
double dot charge stability diagram with current only at 
the triple points, is recovered. 

These different cases represent different coupling 
regimes between the two dots, and between the dots and 
the leads. In the case of V^g = —12 V, the current is 
high whenever the energy level in the left dot is aligned 
with the chemical potential in the left lead. This can 
be understood assuming the coupling between the right 
dot and the right lead is very strong compared to the 
coupling of the left dot to the leads and to the right 




FIG. 3: (a) Measurement of a charge stability diagram around 
one pair of triple points at small bias voltage Vbias = 15 fiV 
and Vbg = -8 V, Vol = Vcr = V, Vcg = -0.5 V and 
T = 120 mK. (b) Schematic drawing of a charge stability 
diagram for two tunnel coupled quantum dots around two 
triple points labelled A and B. The corners of the hexagons, 
at the triple points, are rounded due to tunnel coupling, with 
a deviation from the straight, dotted lines proportional to 
the tunnel coupling strength t. (c) Simulation of a pair of 
triple points for t = 14 /xeV, and the energy independent 
part of the tunneling rates to the leads 7l = 1.26 GHz and 
7r = 1.69 GHz. (d) Maximum current along the cotunneling 
lines in the range of the upper triple point for each value of 
Vgl for the measurement (black crosses) and simulations for 
t = 12 /ieV (green), t = 14 /xeV (blue) and t = 16 /xeV (red). 



dot, and therefore transport is dominated by the left 
dot. For Vqg = —13 V, the opposite situation is real- 
ized, with strong coupling between the left dot and the 
left lead. Fig [2](d), recorded at Vbg = —16 V, corre- 
sponds to a more symmetric situation, where the current 
along the cotunneling lines in both directions is almost 
equally strong. The coupling between both dots Eq is 
very strong compared to the charging energies of the in- 
dividual dots and E^ {E^ ^ 0.5 ■ E^ ^ 0.5 ■ E^), 
leading to almost diagonal lines, which would correspond 
to one large dot delocalized over both islands. Finally, at 
Vbg = —19 V, the situation corresponds again to a well 
defined double dot, with all three tunnel barriers well 
closed {E^ ^ 0.2E^ ^ 0.2E^). 

Fig. [3]^a) shows a measurement for a charge stability 
diagram in the vicinity of one pair of triple points. It has 
been recorded at low temperature T ^ 120 mK and at 
low bias voltage Vbias = 15 /iV to prevent an expansion 
of the triple points to triangles. A negative back gate 
voltage Vbg = — 8 V has been chosen, because only in 
this regime the cotunneling lines are visible even at this 
low bias voltage. A corresponding schematic drawing of 
a charge stability diagram for two tunnel coupled quan- 
tum dots can be seen in Fig. ^h). The tunnel coupling 
leads to rounded edges of the hexagons, with the point 
of charge balance shifted from the original triple point 
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proportional to the strength of the tunnel coupling. Con- 
trary to the situation depicted in Fig. [3l^b), no rounding 
of the corners is visible in the measurement. This allows 
us to estimate an upper bound for the tunnel coupling t < 
20 /ieV. This energy scale is comparable to the tempera- 
ture broadening of the cotunneling lines (/cbT ^ 10 /ieV) 
and about two orders of magnitude smaller than the ca- 
pacitive coupling energy ^1.3 meV. The lever arms 
necessary for the determination of these energy scales 
were extracted from a measurement of the same pair of 
triple points at l^ias = 1 mV. Fig. [3jc) shows the re- 
sult of a numerical calculation of the current based on 
the rate equation using the lever arms and charging en- 
ergies deduced from the measurement. Best agreement 
is found for 7l = 1.26 GHz and 7r = 1.69 GHz for 
the energy-independent part of the tunneling rates to 
the leads, and t = 14 fieV for the tunnel coupHng be- 
tween both dots^^^^^. A more detailed description of the 
calculation can be found in Ref.^^. Measurement and 
simulation are quite similar, except for the finite current 
measured along the line connecting both triple points, 
which is unexpected and can not be reproduced with this 
simple model. Fig. ^d) shows a closer comparison be- 
tween the measurement and the calculation. For each 
value of Vgl in the range of the upper triple point, the 
maximum current in the region of the cotunneling line 
is plotted for the measurement and for calculations with 
t = 12 /ieV, t = 14 /ieV and t = 16 /ieV. The tunnel cou- 
pling determines how fast the current drops as one goes 
away from the triple point along the cotunneling lines. 
Best agreement is found for t = 14 /ieV. However, due to 
the uncertainty in the lever arms in this regime and in 
the electronic temperature, the tunnel coupling strength 
can only be determined up to a factor of two. The differ- 
ence between measurement and calculation in the tails of 
the peak far away from the triple point arises from the 
fact that the peak current from the measurement does 
not take values lower than the noise level. 

Despite the strong capacitive coupling between the 
dots, the tunnel coupling is low. For the observation 
of Coulomb blockade a resistance of the order of h/e'^ is 
required. This resistance may arise from a tunnel bar- 
rier, as it is usually the case for Ga As-based quantum 
dot systems. In graphene, this tunnel coupling may be 
weak if there is a narrow but high barrier separating the 
dots. Such a situation could give rise to strong capacitive 
coupling (see Fig. Efd)) while the tunnel coupling itself 
remains below the experimental resolution. Additional 
resonances in the central constriction^^— and interac- 
tions might lead to an even more complicated situation. 

Fig. [4] displays a closer look at one pair of triple points 
for a finite applied bias voltage of Vbias = i6 mV. This 
measurement is recorded at Vbg = 25 V, a region where 
the dot-lead coupling strengths are weak, and therefore 
no cotunneling lines are visible. Because of the high bias 
voltage, the triple points evolved into triangles^. The 
extent of these triangles gives the lever arms needed to 
determine the energy scales of the double dot system. 
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FIG. 4: Closer look at one pair of triple points at T = 120 mK. 
(a) Current along the detuning axis. The arrows indicate the 
position of excited states, (b), (d): Current as a function of 
Vgr and Vgl for Vbias=6 mV (b) and Vbias=-6 mV (d). The 
dashed line in (b) represents the detuning line chosen for the 
measurement in (a). (c),(e): Corresponding representations 
of the current, numerically derivated by Vgl after smoothing 
over three data points. 



TABLE I: Overview of the main lever arms and energy scales. 
The lever arms and the mutual charging energy Eq are de- 
termined from Fig. [4l^b). The single dot charging energies are 
determined from Fig. [TJa) using the lever arms of this table. 



0.18 



ttGL,L 
0.20 



0.062 



ttGL,R 

0.066 



Eg E^ E^ level spacing 

13.0 meV 14.0 meV 2.4 meV ^ 0.5 meV 



Table [J gives an overview of the main lever arms and 
energy scales. The symmetry of the structure is remark- 
able: the lever arms of both dots, as well as both single 
dot charging energies, are very similar. The charging 
energies are considerably higher than those reported in 
Ref.^^ (single dot charging energies ~ 3 meV, mutual 
coupling energy ^0.4 meV) despite the fact that these 
dots are slightly smaller. This is consistent with the fact 
that in the case of Ref.^^ the sample is partly covered by 
top gates, leading to increased screening. 

Inside the triangles of finite current, additional paral- 
lel lines can be seen. These lines are even clearer when 
plotting the derivative of the current along the VcL-axis, 
taken numerically after smoothening over 3 data points. 
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FIG. 5: dJ/dVoL as a function of Vgr and Vgl for T = 
120 mK and Vbias = —4 mV for different values of in-plane 
magnetic field. c/J/c/Vgl is measured with a lock-in amplifier 
by adding an ac modulation of 200 /iV to Vgl and detecting 
the ac component of the current. 

The most prominent lines run parallel to the baseline of 
the triangles. Along such a line, the detuning between the 
energy levels in both dots is kept fixed. These lines are 
usually attributed to excited states in the right (left) dot 
for positive (negative) bias voltage. The lines can also be 
clearly seen in a cut along the detuning line (arrow in Fig. 
[11(a)). They have a typical level spacing of ~ 0.5 meV for 
excited states in the right dot, and 0.4-0.8 meV in the 
left dot, and are much broader than k^T ^ 10 /ieV as a 
result of inelastic tunneling processes. 

However, at closer inspection additional lines parallel 
to the lower edge of the triangle for ]4ias = 6 mV and 
to the upper edge for Fbias = — 6 mV are visible. Along 
these lines, the alignment between the energy level of 
the left dot and the Fermi energy in the left lead is kept 
constant. In the case of negative bias voltage, these lines 
can not originate from an excited state in the left dot, 
assuming the number of carriers in both dots to stay 
constant. These lines are probably due to modulations 
of the tunneling coupling between the left dot and the left 
lead, because of resonances in this constriction^^. These 
lines are broader than the lines parallel to the baseline, 
and only occur parallel to the nearly horizontal edge of 
the triangle, which corresponds to the direction of the 
stronger cotunneling lines. Parallel to the other edge of 
the triangles, no lines are observable, even when taking 
the derivative in the other direction. 

Fig. [5] presents a study of one pair of triple points for 



different values of the magnetic field, oriented parallel to 
the graphene plane. The figure displays the current dif- 
ferentiated by Vgl, measured directly by applying an ac 
modulation to Vgl, and recording the ac current. Again, 
states parallel to the baseline as well as lines parallel to 
the upper edge of the triangle are visible. The position 
of the triangles in the gate voltage plane almost does not 
change at all up to 5 = 12 T. This is in contrast to the 
case of a perpendicular magnetic field, where the position 
of the triple points and the intensities change significantly 
on a magnetic field scale ^ 250 mT due to the ef- 
fects of the field on the orbital part of the wavefunctions 
(not shown). Two effects of the parallel magnetic field 
on the triple points can be observed. First, with increas- 
ing magnetic field, the number of visible states parallel 
to the baseline increases, which is most pronounced for 
B = 12T [Fig. El^d)]. This effect was observable for 
both pairs of triple points which were studied, and could 
originate from Zeeman splitting. However, it was not 
possible to analyze in detail the appearance of these ad- 
ditional lines because their broadening is similar to their 
spacing. The second effect is the appearance of a line 
parallel to the left edge of the triangle at high magnetic 
fields. This line originates from modulated transmission 
between the right dot and the right lead. The appearance 
of this line with high magnetic field is surprising, as one 
would not expect the in-plane magnetic field to localize 
states. We speculate that some areas of the structure are 
exposed to a finite component of the out-of plane field, 
owing to ripples always present in graphene flakes^!. 



IV. CONCLUSION 

We have studied a graphene double quantum dot in 
different coupling regimes. Despite the strong capaci- 
tive coupling between both dots, the tunnel coupling is 
below the experimental resolution and no roundening of 
the hexagones at the triple points can be resolved. A 
numerical calculation of the current based on the rate 
equation leads to an estimation for the tunnel coupling 
of t ^ 14 /ieV. Inside the finite-bias triangles, additional 
structures can be observed, which we attribute to excited 
dot states, but partly also to imperfections in the tunnel 
barriers. With the application of an in-plane magnetic 
field, additional states become visible within the finite- 
bias triangles. 
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